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Using Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere 
Fluid 

Jeff Schulte, 1 Patrick Kreitzberg, 1 Chris Haglund, 1 and David Roundy 1 

Department of Physics, Oregon State University, Corvallis, OR 97331 

We investigate the value of the correlation function of an inhomogeneous hard-sphere fluid at contact. This 
quantity plays a critical role in Statistical Associating Fluid Theory (SAFT), which is the basis of a number 
of recently developed classical density functionals. We define two averaged values for the correlation function 
at contact, and derive formulas for each of them from the White Bear version of the Fundamental Measure 
Theory functional 1 , using an assumption of thermodynamic consistency. We test these formulas, as well as 
two existing formulas 2 -^ against Monte Carlo simulations, and find excellent agreement between the Monte 
Carlo data and one of our averaged correlation functions. 



I. INTRODUCTION 



There has been considerable recent interest in using 
Statistical Associating Fluid Theory (SAFT) to con- 
struct classical density functionals to describe associating 
fluids 2 11 . This approach has been successful in qualita- 
tively describing the dependence of surface tension on 
temperature. A key input to these functionals is the 
correlation function evaluated at contact, which is re- 
quired for both the chain and association terms in the 
SAFT free energy. The chain term describes the chain 
formation energy in polymeric fluids, while the associ- 
ation term describes the effects of hydrogen bonding, 
both of which can be large. Yu and Wu introduced in 
2002 a functional for the association term of the free en- 
ergy, which included a functional for the contact value 
of the correlation function (described in Section [III E)} 2 . 
which has subsequently been used in the development of 
other SAFT-based functional s 10 ! 11 . Two functionals for 
the chain contribution have recently been introduced, one 
which uses the correlation function of Yu and Wuii and 
another that introduces a new approximation for the con- 
tact value of the correlation function (described in Sec- 
tion IIIID|^ . Given these different approaches, it seems 
valuable to examine this property of the hard-sphere fluid 
through direct simulation, in order to establish the ad- 
vantages and disadvantages of each approach. 



Although these recent works have introduced approx- 
imate functionals for the contact value of the correla- 
tion function 2 -^, there has not been a study that specif- 
ically addresses this contact value for an inhomogeneous 
hard-sphere fluid. In this paper we introduce two defini- 
tions for the locally averaged correlation function of an 
inhomogeneous system. Given these definitions, we will 
present a thermodynamic derivation for each correlation 
function from the free energy functional. We will then 
discuss the correlation functions of Yu and Wu and of 
Gross, and will end by comparing all four approxima- 
tions with Monte-Carlo simulations of the hard-sphere 
fluid at a variety of hard- wall surfaces. 



II. CORRELATION FUNCTION WITH 
INHOMOGENEITY 

We define our terms using the two-particle density 
n( 2 )(ri, r 2 ), which gives the probability per unit volume 
squared of finding one particle at position ri and the 
other at position r 2 . The pair correlation function is de- 
fined by 



#(ri,r 2 ) 
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(ri,r 2 ) 



n(ri)n(r 2 ) 
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In a homogeneous fluid, the pair correlation only depends 
on the distance |i*i — r 2 | and can be expressed as a func- 
tion of a single variable. The contact value of the corre- 
lation function is this correlation function's value when 
evaluated at a distance of the diameter a. It is desir- 
able for reasons of efficiency to limit CDFT functionals 
to one-center convolutions, which leads us to seek a sim- 
plified expression for the contact value of the correlation 
function — which is the same as the contact value of the 
cavity correlation function for hard spheres. In a system 
with an inhomogeneous density, we seek a local value for 
g a . There are two reasonable options for defining such a 
local function: a symmetric formulation (which we refer 
to as S) and an asymmetric formulation (which we refer 
to as A). 

For the symmetric S case, the correlation function at 
contact is given by: 
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where a is the hard sphere diameter and density no is one 
of the fundamental measures of Fundamental Measure 
Theory (FMT), and is ideal for treating touching spheres, 
as illustrated in Figure [TJ 



n (r) = J n(r') 
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This functional gives a value averaged over all spheres 
that touch at the position r. 
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FIG. 1. Set of hard spheres that are included in no(r), which 
consist of those which just touch the point r. 



In contrast, the asymmetrically averaged A correlation 
function is given by 

flftr) = rr~~ tt I ^(r.r + r')^*' (4) 



n(r)n A (r) 
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where the density n^(r) is similar to no, but measures 
the density of spheres that are touching a sphere that is 
located at point r. 



(5) 



Thus g£ corresponds to an average of the two-particle 
density over spheres touching a sphere that is located at 
the position r. 



Fundamental-Measure Theory 

We use the White Bear version of the Fundamental- 
Measure Theory (FMT) functional 1 , which describes the 
excess free energy of a hard-sphere fluid. The White 
Bear functional reduces to the Carnahan-Starling equa- 
tion of state for homogeneous systems. It is written as 
an integral over all space of a local function of a set of 
"fundamental measures" n a (r), each of which is written 
as a one-center convolution of the density. The White 
Bear free energy is thus 



A HS [n] = k B T J ($!(r) + $ 2 (r) + * 3 (r)) dr , 



with integrands 

$i = -n ln(l - n 3 ) 

nin 2 - nyi • n V2 



1 - n 3 
(nl - 3n 2 n y2 ■ n V2 ) 



n 3 + (1 - n 3 ) 2 ln(l - n 3 ) 
367m3 (1 — n 3 ) 2 



(6) 

(7) 

(8) 

I 

(9) 



ra 3 (r) = J 


n(r')e(cr/2 - 


-|r- 


-r'\)dr' 


( ^ / 
n 2 (r) = / 


n(r )o(a/2 — 


1 

|r- 


/ 1 \ 7 / 

r Jar 


n 2 y(r) j 


n(r')5(cr/2 - 


|r- 






ny 2 


n 2 
2ira 


, = — o 

7T(J Z 



using the fundamental measures 
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III. THEORETICAL APPROACHES 

A. Homogeneous limit 

In order to motivate our derivation of the correlation 
function at contact for the inhomogeneous hard-sphere 
fluid, we begin by deriving the well-known formula for g a 
for the homogeneous fluid that comes from the Carnahan- 
Starling free energy. The contact value of the correlation 
function density can be found by using the contact- value 
theorem, which states that the pressure on any hard sur- 
face is determined by the density at contact: 



p = k B Tn contact 
= k B Tng a 



(14) 
(15) 



Since we are interested in the correlation function at the 
surface of the hard spheres, we need to compute the pres- 
sure on that surface. The pressure on a hard sphere 
can be readily computed from the dependence of the 
Carnahan-Starling free energy on hard sphere radius. 



A HS = Nk B T 



4ry - 3r/ 2 



(16) 



where r] = n is the filling fraction. We can thus 
readily compute the derivative of free energy with respect 
to hard-sphere diameter, which is half the force on all the 
hard spheres — since we're changing all their diameters at 
once. To compute the pressure on the spheres, we need to 
divide the total force by the total area of contact, which 
means dividing by NAttct 2 . 
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Using the contact-value theorem, we can thus find the 
well-known correlation function evaluated at contact. 
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Extending this derivation to the inhomogeneous fluid re- 
quires that we find the pressure felt by the surface of 
particular spheres. 
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B. Asymmetrically averaged correlation 
function 

We will begin our derivation of the locally averaged 
correlation function with the asymmetric definition of 
g^(r) given in Equation 2J which is averaged over con- 
tacts in which one of the two spheres is located at posi- 
tion r. This correlation function is related to the contact 
density averaged over the surface of a sphere located at r, 
and can thus be determined by finding the pressure on 
that sphere. We find this pressure from the change in 
free energy resulting from an infinitesimal expansion of 
any spheres located at position r. From this pressure, we 
derive a formula for the correlation function g^(r) as was 
done in the previous section: 



a reasonably smooth density: 

^4/n|-|r-rW)^ (23) 
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a 

In the systems that we study, the density is not reason- 
ably smooth, but we can state empirically making this 
approximation nevertheless improves the predictions of 
our functional, while at the same time reducing its com- 
putational cost by avoiding the need to calculate any ad- 
ditional weighted densities or convolutions. 
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Details of evaluating the functional derivative in Equa- 
tion [21] using FMT are given in the appendix. The equa- 
tion for g£ requires finding convolutions of local deriva- 
tives of the free energy, making this formulation compu- 
tationally somewhat more expensive than the free energy 
itself. 



C. Symmetrically averaged correlation function 

We now address the symmetrically averaged correla- 
tion function, which is defined in Equation [2] This cor- 
responds to the correlation function averaged for spheres 
touching at a given point. In this case, we conceptually 
would like to evaluate the pressure felt by the surface 
of spheres where that surface is located at point r. We 
can approximate this value by assuming that this pres- 
sure will be simply related to the free energy density at 
point r. Through a process similar to the previous deriva- 
tions, this leads to the expression 



1 1 <9$(r) 
no(r) 2 27rt 2 da 



(22) 



where 3>(r) = $i(r) + $2(1*) + $3(r) is the dimension- 
less free energy density. This step is an approximation — 
unlike the analogous Equation [21] — because it assumes 
that we have available a local functional 3>(r) whose 
derivative provides the pressure needed to compute <7<r(r). 
Equation [22] requires that we evaluate the derivatives of 
the fundamental measures n a (r) with respect to diame- 
ter, which leads us to derivatives of the S function, which 
we can simplify and approximate using an assumption of 



D. Gross's asymmetrically averaged correlation 
functional 

One approximation for the correlation function is that 
of Gross^, which is of the asymmetrically averaged variety 
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where ua is the averaged density defined in Equation [5] 
This formula is arrived at by using the density aver- 
aged over all spheres that could be touching a sphere 
at point r in the Carnahan- Starling equation for the cor- 
relation function at contact, given in Equation fT9l 



E. Yu and Wu's symmetrically averaged 
functional 

Yu and Wu developed a functional for the correlation 
function evaluated at contact which is symmetrically av- 
eraged 2 . However, instead of using no as the correspond- 
ing density, they use a density given by 



n Yu (r) = n (r)C(r) 
n 2 • n 2 



(27) 
(28) 



where the function £ is a measure of local inhomogeneity 
at the point of contact, and has the effect of reducing 
this density at interfaces. Because of this difference, the 
correlation function of Yu and Wu cannot be directly 
compared with g^ as defined in Equation [2] Therefore 
in order to make a comparison we move the factors of 
C in Equation [23 from the density into the correlation 
function itself. 
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FIG. 2. Density and correlation function at a hard wall with 
bulk filling fraction of 0.1. 

where gj u is the correlation function as defined in refer- 
ence^, and gj u,s is the function we will examine in this 
paper. 

IV. COMPARISON WITH SIMULATION 

We performed a Monte-Carlo simulation of the hard 
sphere fluid to measure the contact value of the corre- 
lation function for several simple inhomogeneous config- 
urations. For each configuration, we compute the mean 
density, and the contact values of the correlation func- 
tion, averaged as defined in Equations |4] and [2j We com- 
pare these with the functionals presented in sections UlI Bl 
to IIII El We constructed our functionals using both the 
original White Bear functional^ as well as the mark II 
version of the White Bear functional^, but the results 
were essentially indistinguishable on our plots, so we ex- 
clusively show the results due to the original White Bear 
functional. 

We simulate the inhomogeneous hard sphere fluid at 
four hard- wall interfaces. The first and simplest is a 
flat hard wall. We then study two convex hard surfaces. 
One is an excluded sphere with diameter 2<j, which cor- 
responds to a "test particle" simulation with one of a 
hard sphere at the origin with diameter a. The second 
is an excluded sphere with diameter 6<j, which demon- 
strates behavior typical of mildly convex hard surfaces. 
Finally, we study a concave surface given by a hard cav- 
ity in which our fluid is free to move up to a diameter of 
16<j, which demonstrates behavior typical of mildly con- 
cave surfaces. In each case, we performed a low-density 
(filling fraction 0.1) and high-density (filling fraction 0.4) 
simulation. We performed additional computations over 
a wider range of curvatures and densities, but chose these 



FIG. 3. Density and correlation function around a hard sphere 
with diameter the same as those in the fluid, with a bulk filling 
fraction of 0.1. 



as typical examples. 



A. Low density 

We begin by presenting our low-density results, cor- 
responding to a filling fraction of 0.1, which are shown 
in Figures [2j(5l At this low density, the contact value of 
the correlation function in the bulk is only 1.3, indicat- 
ing that correlations are indeed small and that the fluid 
should be relatively easy to model. Indeed, the contact 
density at the hard surface is only around 50% higher 
than the bulk, and the FMT predicted density is close 
to indistinguishable from the true density for each of the 
four configurations. 

The g£ correlation function in each configuration is 
very flat, with only small, smooth changes as the surface 
is approached. Our functional g£ very closely matches 
the Monte Carlo predictions in each case, while that of 
Gross consistently underestimates the correlation at the 
interface by a significant margin. We note that the theo- 
retical curves extend into the region from which the fluid 
is excluded. This value corresponds to the correlation 
function that would be observed in the vanishingly un- 
likely scenario in which there was a sphere present at 
that location. Naturally, we are unable to observe this 
quantity in our Monte Carlo simulations. 

The g% correlation function shows considerably more 
structure, as well as additional variation due to the cur- 
vature of the hard surface. The symmetric correlation 
function is nonzero at locations where spheres may touch, 
which for a convex hard surface means that g% may be 
nonzero in the volume in which hard spheres are ex- 
cluded. In every configuration studied, the agreement 
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FIG. 4. Density and correlation function around a hard sphere FIG. 6. Density and correlation function at a hard wall with 
with an excluded diameter of 6cr, with a bulk filling fraction bulk filling fraction of 0.4. 
of 0.1. 
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FIG. 5. Density and correlation functions near the surface of 
a spherical cavity with diameter 16a at a bulk filling fraction 
of 0.1. 




FIG. 7. Density and correlation function around a hard sphere 
with diameter the same as those in the fluid, with a bulk filling 
fraction of 0.4. 



B. High density 



between the theoretical predictions and the Monte Carlo 
simulation in each case is very poor in the region where 
there should be no contacts at all. Because no is com- 
parable to its bulk value in this region, this means that 
these functionals predict a significant number of contacts 
in the region where there should be none. The correla- 
tion function of Yu and Wu 2 and ours described in Sec- 
tion mEO give similar results, with slightly larger errors 
in our prediction. 



At a higher density corresponding to a filling fraction 
of 0.4, correlations are much stronger, with the bulk con- 
tact value of the correlation function of 3.7, as seen in 
Figures [6E1 This results in larger oscillations in the den- 
sity at the hard surfaces, and correspondingly more in- 
teresting behavior in the correlation function near the 
interface. The density predicted by the White Bear func- 
tional agrees reasonably well with the simulation results, 
although not so well as it did at lower density. The dis- 
crepancies are largest in the case of the spherical cavity, 
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FIG. 8. Density and correlation function around a hard sphere 
with an excluded diameter of 6cr, with a bulk filling fraction 
of 0.4. 



FIG. 9. Density and correlation functions near the surface of 
a spherical cavity with diameter 16a at a bulk filling fraction 
of 0.4. 



in which the DFT considerably underestimates the range 
of the density oscillations. 

The asymmetric version of the correlation function 
once again displays relatively smooth behavior with a 
few small oscillations near the interface, and a somewhat 
elevated value within a diameter of the hard surface, with 
the magnitude of this elevation somewhat different in 
each configuration. As was the case at low density, our 
correlation function g£ matches very closely the Monte 
Carlo data, reproducing quite well the structure near the 
interface in each configuration, although in the spher- 
ical cavity (Figure [9j), there is a small, but significant 
discrepancy, comparable to the discrepancy found in the 
density itself. In each case, the correlation of Gross dra- 
matically underestimates the correlation at the interface, 
at one extreme by 40% in the case of the spherical cav- 
ity (Figure |5J), and at the other extreme by 15% in the 
test-particle scenario (Figure [7]). 

The symmetrically averaged correlation function shows 
considerably more structure near the interface at high 
density, and this structure varies considerably depending 
on the curvature of the hard surface. In each case, this 
structure is not reflected in the theoretical predictions, 
neither that of this paper, nor that of Yu and WiA As 
was the case at low density, both functionals give signif- 
icant and finite values in the region in which there are 
no contacts, but at high density they also miss the large 
oscillations that are present near the fiat wall and the 
concave surface (Figures [6] and [9j). As was the case at 
low density, the functional of Yu and Wu^ gives slightly 
better agreement with the simulation results than that 
which we derive in Section IlII CI 



V. CONCLUSION 

We investigated several approximations to the contact 
value of the correlation function for inhomogeneous fluid 
distributions corresponding to flat, concave, and convex 
walls. We defined and simulated two averages of the cor- 
relation function, an asymmetric A average centered at 
the location of one of the two spheres that is in contact, 
and a symmetric S average centered at the point of con- 
tact of touching spheres. For each average, we derived a 
functional form from FMT, and also found an approxi- 
mation that has been used in the literature. When com- 
pared with essentially exact Monte Carlo simulations, the 
A correlation function derived from Fundamental Mea- 
sure Theory in Section IIII Bl gives excellent results for 
each surface, at both high density and low density. The 
other three approximations that we studied all showed 
significant and systematic deviations under some circum- 
stances. Thus, we recommend that creators of SAFT- 
based classical density functionals consider using the g£ 
functional defined in Section HlIBi 



APPENDIX 



The expression for the asymmetric correlation function 
<7^ (r) (Equation |2T]) involves the functional derivative 
f^y • In this appendix we will explain how this deriva- 
tive is evaluated. We begin by applying the chain rule in 



the following way: 
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This expression requires us to evaluate anc ^ n ° recm ce to sums of the terms above: 



. The former is straightforward, given Equa- 
tions [7]l9l and we will write no more about it. The func- 
tional derivatives of the fundamental measures, however, 
require a bit more subtlety, and we will address them 
here. 

We begin with the derivative of ns, the filling fraction, 
which we will discuss in somewhat more detail than the 
remainder, which are similar in nature. Because the di- 
ameter <j(r) is the diameter of a sphere at position r, we 
write the fundamental measure 77,3 (V) as 



n 3 (r') = |n(r")6 (^p - |r' - r"|) dr" 



(32) 



where we note that <r(r") and n(r") are the diameter 
and density, respectively, of spheres centered at posi- 
tion r". Thus the derivative with respect to the diameter 
of spheres at position r is 



Sn 3 (r') = 1 
Sa(r) 2 
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The vector-weighted densities nyi and ny2 give terms 
analogous to those of rt\ and 



invi(r') _ »W. yWr)/2 _ |r ,_,. r r-r' 
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Thus there are four convolution kernels used in comput- 
ing one scalar and one vector delta function, and one 
scalar and one vector derivative of the delta function. 



This pattern will hold for each fundamental measure: 
because we are seeking the change in free energy when 
spheres at point r are expanded, the integral over den- 
sity is eliminated. To compute the correlation funtion 
g£, we convolve this delta function with the product of 
the density and a local derivative of <£(r): 



SA HS f ggKrQ ( , w /9 , , |W , 



(35) 



As we shall see, there are only four convolution kernels, 
leading to four additional convolutions beyond those re- 
quired for FMT. 

The functional derivative of ri2 introduces our second 
convolution kernel, which is a derivative of the delta func- 
tion. 



6n 2 (r') _ 1 
Sa(r) ~ 2 



^\a(v)l2-\v' -v\) (36) 



The derivatives of the remaining scalar densities rt\ and 
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